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Abstract 

We compute the deeply virtual Compton scattering (DVCS) am- 
plitude for forward and backward scattering in the asymptotic limit. 
Since this calculation does not assume ordering of the transverse mo- 
menta, it includes important logarithmic contributions that are be- 
yond those summed by the DGLAP evolution. These contributions 
lead to a power-like behavior for the forward DVCS amplitude. 
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1 Introduction 



The deeply virtual Compton scattering (DVCS) has received attention re- 
cently since it provides a new and interesting theoretical framework for per- 
turbative and nonperturbative QCD applications. |], |2|, ^] The basic mecha- 
nism of this process is reminiscent of the usual deeply inelastic (DIS) lepton- 
hadron scattering process which begins with the absorption of a virtual pho- 
ton by the quark constitutient. The DVCS differs from DIS in the second 
step when the same quark constitutient then emits an on-shell photon. The 
kinematics of the DVCS reaction are characterized by a high invariant energy 
and a large virtuality of the initial photon compared to the quark mass. 

The DVCS amplitude is, in a sense, a generalization of the conventional 
deep inelastic structure (DIS) functions, and can be expressed using asym- 
metric distribution functions. To incorporate the radiative corrections into 
the DVCS amplitude, the conventional DGLAP evolution equations @] have 
been used 0, [3[] (for review on the DGLAP evolution for the DVCS see || 
and references therein). As a typical for the DGLAP result, the exponential 
asymptotical high energy behavior of the DVCS amplitude was obtained. 
Such asymptotics are true for the kinematical region where the invariant to- 
tal energy and the virtuality of the off-shell photon are of the same order, i.e. 
in so-called the "hard" kinematics. However, in order to obtain the DVCS 
asymptotics in the other kinematical regions, and in particular in the "semi- 
hard" (or the Regge type) one, with the invariant total energy being much 
greater than the off-shell photon virtuality, a generalized set of evolution 
equations is required. 

Although predictions based on the DGLAP evolution yield good agree- 
ment with present experimental deep inelastic scattering data,|6|, 0, |[ from 
the theoretical point of view this approach is insufficient when the kinemat- 
ics are of the Regge type where the total invariant energy is much greater 
than the other invariants, including the virtuality of the deeply virtual pho- 
ton. In this kinematic region, the DGLAP equations can not account for 
contributions independent of this virtuality. Taking these contributions into 
account for the DIS structure functions led to a power-like dependence on 
the total energy both for photon scattering from an unpolarized quark (the 



BFKL PomeronH) and for the scattering from a polarized quark. [10 



Naturally, one can expect a similar change in the asymptotic behavior of 
the DVCS amplitude in the Regge limit. To demonstrate this behavior, let us 
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Figure 1: The fundamental DVCS process: scattering of a deeply virtual 
photon (qi) off a nearly on-shell light quark (pi) to produce an on-shell photon 
(q 2 ) and a quark (p 2 ). 

consider the scattering of a deeply virtual photon off a nearly on-shell light 
quark which produces a photon and a quark — both on-shell, (c/., Fig.l). In 
Fig.l, s , t and u are the standard Mandelstam variables: 

s = (pi + qi) 2 

t = fe-gi) 2 (l) 

u = (q2~Pi) 2 ■ 

In the Born approximation, only the two Feynman graphs in Fig. 2 con- 
tribute to the DVCS process: Beyond the Born approximation the value 
of the QCD-radiative corrections strongly depend on the kinematics of the 
scattering. In particular, in the two regions 

s ~ — u 3> —t Forward Region (2) 

s ~ —t^> —u Backward Region (3) 

the QCD radiative corrections are large, whereas in the region s ~ — t ~ — u 
they are comparatively small. In the center of mass system (CMS), s ~ 
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Figure 2: Born graphs for the DVCS process. 

—u ^> —t corresponds to the Forward Region (Eq. (^j)) (small scattering 
angles) , and s ~ — t ^> — u corresponds to the Backward Region (Eq. 
Both of the above regions are of the Regge type. When 

s » Q\ (4) 

where Q 2 = —q\ is the photon virtuality, the Regge theory predicts a power- 
like s-dependence for the DVCS amplitude in each of the regions. Such 
behavior is result of summing the QCD-radiative corrections to all orders in 
the QCD-coupling. 

Below we obtain the asymptotical behavior of the DVCS amplitudes for 
the "semi-hard" kinematical regions (Eq. (0),Eq. (f|); Eq. ([|),Eq. (^)) in the 
leading logarithmic approximation (LLA) by solving the Infrared Evolution 
Equations (IREE) |TT|, [HJ for them. As in the case of the DGLAP evolution, 
in order to obtain the evolution in the Regge regions above, we consider the 
DVCS from an on-shell quark. Recently, the IREE method has been used 
to obtain the small-x asymptotic limit of the DIS spin-dependent structure 
functions [10] . While references jl], [| [| include consideration of the hadronic 



component, we study only the partonic part of the DVCS which can be 
investigated with the perturbative methods of QCD. Since we consider the 
DVCS from a quark, our results depend on the factorization scale \i\ this 
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factorization scale dependence will cancel out in the expression for the DVCS 
from a hadron. Operationally, within the LLA-accuracy, this effects the 
replacement of /i by some mass factor of order of the hadron mass. 

The paper is organized as follows: In Sect. 2 we briefly review the IREE 
method we shall use to obtain the expression for the DVCS amplitude in 
the forward region, Eq. (|j), Eq. (|]). In Sect. 3 we investigate the DVCS 
amplitude in the backward region, Eq. (Q), Eq. ([|). In Sect. 4 we calculate the 
non- forward structure functions for DVCS in the Regge region Eq. (||), and 
express them in terms of the Reggeons corresponding to on-shell scattering. 
Finally, in Sect. 5 we discuss the implications of these results. 



2 Forward DVCS 

Our calculation will follow the general prescription of the IREE method given 



in Ref. [O, 12 1. All the virtual momenta, ki, can be expressed through the 



Sudakov parameterization as follows: 

ki = otipx + (3iq[ + k ± , (5) 

q{ = 9i - 7r~^ — Pi » ( 6 ) 

2j?i • qi 

k±-pi = k ± • qx = , (7) 

where we have neglected the masses of incoming and outgoing quarks. 

In order to avoid infrared singularities, it is necessary to introduce an 
infrared cut-off /i in transverse momentum space so that 

ki > pi . (8) 

In the LLA, such regularization does not violate gauge invariance; there- 
fore the Gribov bremsstrahlung theoremfO is applicable in this case. (Cf., 



Ref. [fl4| , |lq|.) The essence of this theorem is that the virtual gluon with the 
minimal k± can be factorized out of the amplitude. Consequently, the prin- 
ciple contribution comes from graphs where the gluon propagator is attached 
only to external lines. Additionally, the integrations over the remaining glu- 
ons do not depend on the longitudinal momentum of this factorized gluon. 
There is a dependence on the factorized gluon's transverse momentum — this 
k± becomes a new infrared cut-off for the remaining integrations. 
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Since Q CD-radiative corrections in the LLA do not violate the spin struc- 
ture of the Born approximation, the forward DVCS amplitude, Mp, is given 
by: 

Mp = M|°™ $ 

where we define 

Z\ = ln(s//i 2 ) , 
z 2 = ln(Q 2 //, 2 ) , 

z 3 = ln(-t//i 2 ) . (10) 

Although the value of fi 2 is not fixed, we presume it is much less than 
any of s, Q 2 , or —t. Instead of setting fi 2 to a fixed value, let us differentiate 
Mp with respect to fi to obtain our IREE: 

2 dM F dM F dM F dM F . . 

The LHS of Eq. (0) provides half of the IREE we desire for Mp, (c/., Fig.3). 
In order to obtain the remaining half of the IREE, we should look for the 
virtual quark or a gluon with the minimal transverse momentum. 

If the particle with minimal transverse momentum is a gluon, then it can 
be factorized, (graph (b) in RHS of Fig.3). Integration over the longitudinal 
momentum of the gluon yields (after differentiation over ln(/x 2 )), 

-2\z 3 M F (12) 

with 

A = a s C F /A7i , (13) 

C F = (N 2 - 1)/2N, and N = 3. 

If the particle with minimal transverse momentum is a quark, then the 
process is more involved. In order to evaluate the leading-log (LL) contribu- 
tion, there should be another quark on the ladder with minimal k± from the 
same ladder rung (graph (c) in RHS of Fig.3). However, as long as we keep 
—t^> fi 2 , the LL contribution arises only when k\ > —t. Then the integra- 
tion region over k± does not include \x as the lower limit when —t^> fi 2 , and 



ln( S //i 2 ),ln(Q 2 /^ 2 ),ln(-t//i 2 
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Figure 3: The general form of the evolution equations for DVCS. 



the whole term (c) in Fig. 3 will be independent of fi, and any derivatives 
w.r.t. will vanish. 

By the same reasoning, we can dismiss the case where, instead of quarks, 
two ladder gluons have minimal values of k±, (corresponding to graph (d) in 
Fig.3). 

The Born contribution (graph (a) in Fig.3) also does not depend on \i 
when —t ^> fi 2 , and this term also vanishes when differentiated w.r.t. /i. 

We therefore find the RHS of the IREE is given solely by graph (b) in 
Fig.3, and we obtain 



dMp dMp dMp % , 
+ + = -2\z 3 M F 



dz\ dz2 dz 3 
The solution of this IREE is: 



(14) 



{z x - z 3 , z 2 - z 3 ) e Az 3 



(15) 
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with to be defined below in Eq. (|T7D. 

We will discuss now the kinematic region where 



s > Q 2 > -t . (16) 

We then have —t = fi 2 , which implies z 3 = ln(— t/fi 2 ) = 0, and consequently 
4>(zi — z 3 , z<i — z 3 ) — > <p(z\, z<i). In this kinematic limit we define 

Mf = <f>(zi,z 2 ) (17) 

where Mp is the DVCS amplitude for the scattering at (nearly) zero angle, 
i.e. for —t = fi 2 . We define the function using the matching condition of 
Eq. (|TTD; therefore, we require an independent calculation of Mp to obtain 
0. This amplitude does not depend on z 3 . Thus, we have re-expressed the 
scattering amplitude in terms of a simpler set of variables, and eliminated one 
variable (namely, the smallest one). Let us note that when t ^ 0, the DVCS, 
even when regarded as the elastic 2^2- process, is always accompanied 
by bremsstrahlung consisting of photons and gluons with energies less than 
yf—i. A part of them have energies less than fi. Such very soft radiation 
does not affect kinematics of the DVCS. In the LLA it would factorize out of 



the scattering amplitude (for details see |l4j]), and would cancel with virtual 
gluon/photon contributions of the same energy scale in expressions for the 
cross sections. Through the present paper we neglect both the emission 
and the radiative corrections for bremsstrahlung with energies below [l. So, 
within LLA-accuracy we identify Mp as the usual DVCS-amplitude with 
t = 0. While we begin our calculation with a complicated process involving 
off-shell exchanges, our final result involves elastic scattering amplitudes with 
only on-shell particles. For the DVCS process, the relevant process is the 
elastic quark-photon scattering amplitude with all the particles on-shell. 

We now construct the IREE for Mp as defined in Eq. (|17|). Recall that 
we are working in the limit where z 3 = ln(— t/fi 2 ) = 0. In the following, we 
will show that the only contributions to Mp come from graphs (a), (c), (d) 
of Fig. 3. Let us now demonstrate this result. When the incoming and the 
outgoing quarks in Fig.l are polarized, there are spin- dependent, as well as 
spin-independent, contributions to the DVCS amplitude. In the Born ap- 
proximation, the spin-independent contribution is symmetric with respect to 
permutations of s and u, whereas the spin-dependent contribution is anti- 
symmetric. This is also true when radiative corrections in the LLA are taken 
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into account. Therefore, the spin-independent part of the DVCS amplitude 
has a positive signature, the spin- dependent part has a negative signature. 

Consequently, we find it convenient to decompose the DVCS amplitude 
into a spin-dependent amplitude N F , and a spin-independent amplitude Up: 

Mp = Np + Up (18) 

We make the same decomposition for the Mp, the DVCS amplitude in the 
z 3 = limit: 

Mp = Np + Up (19) 



The t-dependence of Mp is given by Eq. (|T5| ) provided Mp is known. 

Let us first discuss the spin- dependent part, Np, of the DVCS amplitude 
Mp. If we neglect mass of the initial/final state quark, the quark spin is 
then collinear to its momentum. Since the transverse spin contribution is 
proportional to the quark mass, it is negligible in this limit. 

For simplicity, we consider the case where the Q CD-radiative corrections 
to the Born approximation correspond to adding extra gluon propagators 
to both graphs of Fig. 2. In this case, term (d) of Fig.3 (with a two-gluon 
intermediate state) does not contribute to to the RHS of IREE. Neither terms 
(a) and (b) also do not contribute because t = 0. 

We implicitly define F via a Mellin transform 

rtco A, ) _ 

N F {z u z 2 ) = / i{uo)F{uo,z 2 ) (20) 

J —too ZTTt 

with the negative signature factor 

, . e~ l7ruj - 1 -ma; . . 

£M = — 2 — ' (21) 

The IREE for the spin-dependent part takes the following form: 

(A. + JjF(uj,z 2 ) = ±F(u,Z2)f (u) (22) 

where fo is the Mellin amplitude for quark-quark elastic scattering with all 
quarks on-shell. It corresponds to the lowest blob of graph (c) in Fig.3., and 



explicit expression was obtained in Ref. JTT 



.1/2 

fo = An 2 {u 



2vr 2 V 2-k 2 uj 



[23) 
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with 

^(w) = g 2 N ^(e^IK^y)) , (24) 

and 

y = , (25) 

where g = \/47ra^ is the QCD-coupling (which is fixed in the LLA), and 
D -i is the parabolic cylinder (or Webber) function. 

2N' 1 

The solution of Eq. (|2"2|) is 

za) = e-^-/"/*^)^ . (26) 



Thus, 



toe 



, / x \ ( s \ (Q 
Eq. (|27D fixes the anomalous dimension 



/o(")/8*r 



2 



= / ^MEHb 1^1 (27) 



tot 



7 = / (^)/87T 2 (28) 

for iVp. The unknown function <p{uj) can be fixed by the matching condition: 

Nf = C{z x ) (29) 

when Z2 = 0. The new amplitude C{z\) in Eq. fl2T|) is the amplitude of 
the elastic Compton scattering with both photons on-shell. The IREE for 
C(zi) is obtained by canceling the ^-dependence in the LHS, and adding 
to RHS the contribution of term (a) in Fig. 3 which now is /^-dependent. 
Expressions for the anomalous dimensions of the spin-dependent structure 
function g\ incorporating all nonleading order contributions were obtained 
in Ref. |l2j in the double logarithmic approximation. The resummation of 
the 0(a™ +l ln 2n x) term in the singlet evolution equation was performed in 
Ref. 0. 

Combining the above formulas Eq. (fOp, Eq. (|T7D, Eq. ([TJ)-Eq. ([2"5]), we 
arrive at the result for the spin- dependent part of the DVCS amplitude: 

g 2 C F J^oo 2m\Q 2 ) lu - f /8ir 2 \^ ) 

(30) 
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where e 2 = y/Ana. To identify the contribution to the spin-dependent DVCS 
amplitude from the Born process (term (a) in Fig. 3), we neglect the difference 
between the quark mass and \i: 

sf- = r %l (ix t m (3D 

J -too 2m \fi 1 J uj 

The asymptotic high energy behavior in the forward region, Eq. (|2p, for 
N F is (c/., Ref. 0): 




(32) 



with 



1 + 2A^ ■ (33) 



7T 



We can incorporate term (d) of Fig. 3 into the IREE for Fn as was done 
in in Ref. [K|. The result is to change a in Eq. ( |33"D to 



' - H 2 ~V T (34) 

We now discuss Uf, the spin-independent part of Mp- Thus, we are inter- 
ested in the portion of the amplitude with the positive signature. For t = 0, 
term (b) in Fig.3 does not contribute. In contrast to the spin-dependent part 
of the DVCS amplitude, the contribution of term (c) is now small compared 
to term (d) by a power of s; therefore we neglect term (c). The negligible 
contribution of term (c) can be anticipated since the Regge theory predicts 
the asymptotic behavior of the scattering amplitude to be ~ s (ii+i2+-)- n + 1 ) 
where n is the number of t-channel intermediate particles exchanged, and j'j 
are their spins. 

Repeating the spin-independent calculation as in the previous spin-dependent 
case, we obtain the asymptotic behavior of Up in the Regge limit of the for- 
ward region, Eq. @): 

Up ~ i?( S ,Q 2 )e-^ ln2 (-^ 2 ) (35) 
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where R(s, Q 2 ) is the amplitude for Compton scattering. Since this process 
is mediated by the Pomeron, we obtain 

where Ap is the Pomeron intercept, and jp is the Pomeron anomalous di- 
mension. 

Writing the virtual photon momentum as q± = q' — (p DVCS amplitude 
can be expressed in terms of the asymmetry parameter (, Q 2 = (s. 

3 Backward DVCS 

The backward DVCS process is straightforward to analyze given the previous 
forward DVCS calculation. Since, in the Born approximation, the contribu- 
tion of graph (b) is much greater than graph (a), we can neglect graph (a). 
Furthermore, neither graph (c) nor graph (d) contributes to the IREE in the 
backward region. 

Recall that the kinematic configuration in the backward region is specified 
by (Eq. (D): 

s ~ —t^> — u Backward Region (36) 

Therefore, intermediate partons with minimal value of k± can only exist 
in the w-channel (when u ~ 0). The -u-channel two particle intermediate 
state can not consist of two quarks or two gluons due to fermion charge 
conservation. Instead, a new contribution given by graph in Fig.4 can now 
contribute. However, it is known that this contribution is beyond the LLA. 
Therefore, only terms (a) and (b) of Fig. 3 contribute to the RHS of the IREE 
in the backward region. 

Finally, we note that there is no difference between asymptotic behavior 
of the polarized and unpolarized parts of the backward DVCS amplitude Mb- 
M B depends on the variables z 1 , z 2 , and a new variable z 4 defined as: 

z 4 = ln(-u//i 2 ) . (37) 

For backward scattering z\ w z 3 , so these variables are not independent. In a 
similar manner, there was no ^-dependence for the forward DVCS amplitude 
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Figure 4: The two-particle intermediate state for the backward DVCS pro- 
cess. 

because z± w z\ in this region. The IREE for Mb is (c/., Eq. (|14|)): 
8Mb dM B dM B 

-7, h -5 h -7, — = -2Xzi . (38) 

OZ\ OZ 2 OZi 

with A defined in Eq. ([Hf). The solution for Mb is 

M B = Mi orn $ B (zi-z 4 , z 2 -z A ) e~ Aln2 21 . (39) 

If we assume the hierarchy s > Q 2 > u, we can fix the unknown function 
with the matching condition: 

M b = Mb for 24 = , (40) 
12 



where Mb is the DVCS amplitude for the backward scattering of a collinear 
incoming quark and outgoing gluon. It satisfies 

which yields the solution 

Mb = $ B { Zl -z 2 ) e~ xln2zi (42) 

Finally, we can specify When z 2 = {Q 2 = /i 2 ), M B must coincide 
with the amplitude Cb for the backward Compton scattering. Cb was cal- 



culated in the LLA by direct calculation of the relevant Feynman graphs. JIT) 
The IREE for C B is: 

^ = -2\ Zl C B (43) 

with the boundary condition 

C B = Cf °™ for z x = (s = /i 2 ) . (44) 
Combining the above results, we obtain at the final expression: 

M B = M Born e -^C F ln\sm ^ 



4 Non-forward structure functions 

The DVCS structure functions are a generalization of the conventional DIS 
structure functions. The novel feature DVCS structure functions is the ex- 
plicit dependence on the fraction ( of the longitudinal momentum r trans- 
ferred to a hadron, r = (p. The Q 2 evolution of these structure functions is 
treated in the LLA in the framework of a DGLAP type of equation^]. Here 
we consider the non-forward quark structure function in the DLA, which we 
have applied above to the forward and backward DVCS scattering ampli- 
tudes. 

We can factorize the amplitude of the DVCS process and express the 
quark structure function as a four-quark amplitude Mij tSp (p, k, k — r), (cf., 
Fig. 5). The two lower quarks in the amplitude Mij jSp carrying the momenta 
p and p — r are on the mass shell, while two upper ones with the momenta k 
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and k — r are virtual. The spinor indices {i,j} belong to the lower lines, and 
they are summed with the indices of the initial and final quark spinors Uj (p) 
and Ui(p — r). The total amplitude is given by the integral of the function: 

M sp (p, k,k-r) = Ui(p-r) M ijiSp (p, k,k-r) Uj(p) 

with the elementary quark-photon scattering amplitude: 

d A k 



rpfiu 



(2*) 



Tr M(p, k,k-r) x 



x 



k 



k + q 



fc — r fc „ . 

1 T — 7^ 



k — q' 



k — r 



(k + q) 2 1 (k-r) 2 k 2 1 (k-q') 2 1 (k - r) 2 
The trace in T^ u is taken over the upper line indices {s,p}. 




Figure 5: The graphs for the non-forward structure functions 

To derive the IREE for the amplitude Mij !Sp (p, k,k — r), we decompose 
this using a complete set of 16 gamma matrices: 

16 

Mi j>sp = E 7,1 1% Mab- 

A,B=1 
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where the sum runs over the basis set 7* = {1, 7 A , c Xcr , 7 A 7 5 , 7 5 }- 

In the double logarithmic approximation (DLA), the scalar function Mab 
is comprised of a product of the Born terms which can depend on the invariant 
energy s = (k — p) 2 and the virtualities l\ = k 2 , l\ = (k — r) 2 . Note the 
u — (k + p — r) 2 invariant is not independent for t ~ since s + u = l\ + l\. 

The essential contribution of the amplitude M j; - SJ , comes only in the region 
where both virtualities {/ 2 ,/f} are negative and of the same order, i.e., l\ ~ 
l\ ~ I 2 = — 1\. Therefore, it is natural to write il% sp as a sum of an s-channel 
and a u-channel piece: 

A,B V" " / 

+ \ Y.(r%,(l B )j, Fab (tT'S) ' 

A,B \ P P / 

It is easy to verify that there is no double logarithmic (DL) contribution 
to the IREE in the kinematic region I 2 7^ since Z 2 plays the role of the 
infrared cut-off (instead of fi 2 ). Thus, 



d 8 
ds dz 



s^r + ^\ F S AB = 



where we define 



. i 2 

z = In — 

fi 2 



In terms of a Mellin transform, 

F AB = I^A-A f b AB^,z) 




with a corresponding definition for -u-channel function F AB . These functions 
satisfy the differential equation: 

d rS,U _ rS,U 

~q~ z jab - v Jab ■ 



15 



We can write the solution as: 



■UJZ 



S U 

where the function R^ B is determined by the boundary conditions in the limit 
where virtualities become vanishingly small, I 2 ~ fi 2 which implies z — 0. 
This kinematic region is the Regge limit which describes the scattering of 
on-shell quarks with approximately zero momentum transfer. 

The symbol R A ' B is the coefficient of the Regge amplitude when expanded 
in the {7^,7^} basis. We can introduce an amplitude ^R^g{ which cor- 
responds to the amplitude Mj 3 - jSp in the Regge limit, and expand this in the 
{7 , 7 B } basis as follows: 



A S££ r = Y.u%')ur{p)T^u i ] {p)ul{ q ') 

ijsp 

= - E U h (q'h A U h '(p) ■ U f '(ph B Uf(q>) R S AB 
S A,B 

- E U f (q'h A U h '(p)-U f '(ph B U h (q') R% , 



+ 



U A,B 



where h, f, h', and /' are the quark helicities. 

The spinor structure of the Regge amplitude is determined by the Born 
terms, so we have 

^ijjsp 7si7jp 9 Rs ~t~ 7pi7js Ru 

s u 

In the Regge theory framework, it is more natural to work with amplitudes 

of definite signature rather than s-channel and u-channel amplitudes, rather 

than with the forward and backward ones, therefore, we work with i? 1 * 1 defined 

by: 

R s = R(+) + R(~) 

Ru = i?« - R^ . (47) 

The form of the Regge amplitude dictates the final structure for the four- 
quark amplitude: 
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To extract the non-forward structure function, the four-quark amplitude 
must be substituted into the total amplitude. For t ~ 0, the total amplitude 
is given by the sum of the two transverse tensors, 



2p ■ q 



7 9T T, 



2p ■ q' 



IE 



To 



p ■ q' 



where T\ is symmetric, and T2 is antisymmetric. We can express T\ and T2 
in terms of the quark (q) and antiquark (q) non-forward structure functions 



Jf q {x) and Gl' q (x) 



Ti = U{p — r)q'U{p) 
T 2 = U{p-r)q'j 5 U{p) 



dx 



dx 



x — ( + i5 
1 



+ 



x — i5 
1 



x — ( + i5 



x — i5 



(Gj(x) + G|(x)048) 



(Here we follow the definitions given in Refs.Q.) 

Introducing the Sudakov variables, the explicit expressions for T\i are: 



with 



T-i 



To 



d A k „ „ 



x 



(2?r) 4 « 
1 



Tr M(p, k,k — r)kq'(k — f) 



1 



k 2 (k — r) 



+ 



l3-( + i5 /3-i5 
d 4 k „ „ 



1 



x 



1 1 



P - C + iS P-iS 



d A k 



— da dp dkj 
2p-q' . 



(49) 



The integrals over (3 are restricted to the interval < P < 1 which follows 
from the position of the poles in the a-plane. Examining these expressions 
demonstrate that to find the structure functions, the amplitude M must be 
integrated over the virtual momentum k while keeping (3 fixed. In performing 
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this integration, we need only keep the terms of the form da /a dk\jk\ for 
DL accuracy. These terms come from the region where 

k 2 = a(3s' - k\ « - k\ 

{k-rf = a(P-C)s' - k\ « - k 2 ± 

s = (k-p) 2 = a{(3-l)s - k]_ « a{(3-l)s 

u = (k+p-r) 2 = a(P + l-t)s' a(p + 1 - CK- 

These relations simplify the process of extracting the Regge limit. First, 
they imply that both virtualities {I 2 , 1 2 } in the four-quark amplitude must 
be of the same order, i.e., l\ ~ l\ ~ — this justifies the assumption that 
was adopted previously. Second, only the quadratic terms in k±_ from the 
spinor trace are necessary for the non-forward amplitudes T lj2 . We can then 
compute these traces as follows: 

Tr M k ± q' k ± = U(p - r) 7 A k ± q'k ±1 B U(p) {-F S AB + —F% B ) 

s u 

= -2(-R s + -Ru) k 2 ± U(p - r) q' U(p) 
s u 

Tr M k^ 5 k ± = U(p - r) 7 A k ± q'l 5 k±l B U(p) {-F S AB + -F AB ) 

s u 

= 2(-R s + -Ru) k 2 ± U(p - r) g'7 5 U(p) . 
s u 

To calculate the structure functions, we must replace the infrared cut-off 
/i 2 in the Regge amplitude by — kj_, and then integrate over the loop mo- 
mentum. The forward and backward amplitudes will yield the quark and 
antiquark functions, respectively. In the DL region, the transverse momen- 
tum is 

k\ > Max{\ap\s' , \a(P - ()\s'} . 

The remaining integral over a is purely logarithmic, and the symmetric in- 
tegration limit yield in. 

Finally, working with DL accuracy we retain only the non- vanishing terms 
in the limit uj — > to arrive at the following structure functions: 



(50) 
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8tt 5 



C h 



1 + /5-C 



duj . . 1 
— ifr(w) - x 
Z7T u; 



' l+/3-C ^ 

/5 , 



(51) 



The 



The explicit expressions for R$, Ru can be taken from Refs. [XT], 
functions are expressed by Regge amplitudes using the same formula. 

The integral over u runs along imaginary axes to the right of the singularities 
of the Rs,u functions. Non-zero contributions exist only when the arguments 
in the brackets are larger than unity, so that the actual support of the first 
term on the RHS of Eq. Q50|) and Eq. Q5Tp is C/2 < (3 < 1/2. The Regge 
limit implies the argument of the Mellin transform is be asymptotically large, 
which corresponds to the the region where £ <C 1 and /3 ~ £• 

The absence of a Q 2 dependence in these expressions resembles the DGLAP 
type of evolution equation for the non-forward structure functions in Q 2 — > oo 
limit. In contrast to the DGLAP evolution equations, the DVCS evolution 
equations are not restricted to the region (3 < (. 

Notice that for j3 = 0, we find that JF' ,,? (0) ^ in contrast to the result 
for the DGLAP non-forward functions. The Mellin moments M% ,q (Q = 
Jq dx x n ~ x \x) for the integer n are not necessarily polynomials of degree 
n with respect to the transfer momentum fraction (, as in the DGLAP case. 



5 Conclusion 

We obtained the explicit expressions for the DVCS amplitudes in the Regge 
limits of forward and backward scattering. For the forward scattering region, 
the asymptotic amplitude takes the form of the DIS structure functions times 
a Sudakov exponential factor of the form exp[— (a s /4ir) Cf ln 2 (— t/fi 2 )]. The 
DVCS amplitude in the backward scattering region is purely of the Sudakov 
type with an exponential factor of the form exp[—(a s / An) C^ln 2 (s/yLi 2 )]. We 
have not used the DGLAP evolution equations which assume the transverse 
momenta to be strictly ordered in virtuality. This enables us to take into 
account important logarithmic contributions that are beyond of the realm 
of the DGLAP approach. Specifically, there is no ordering of the transverse 
momenta in the Regge region, and this allows the internal transverse mo- 
menta to exceed Q 2 , in contrast to the DGLAP method. Such contributions 
lead to the power-like behavior of the DVCS amplitude. 
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